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A theoretical model is presented to explain the spin gap 
observed for CaV40g. The underlying lattice of the 1/5- 
depleted square lattice favors a formation of plaquette res- 
onating valence bond state. Inclusion of the frustrating sec- 
ond neighbor interaction enhances this tendency, leading to 
a quantum disordered state of a two dimensional spin-1/2 
Heisenberg model with a sufficiently big spin gap compatible 
with experiments. 

PACS numbers: 75.10.Jm, 75.30-m, 75.40-s 



Quantum disordered phases with a spin gap are of 
great current interest. This topic has gained additional 
momentum by Anderson's proposal of the resonating va- 
lence bond (RVB) state in the undoped parent materials 
of high temperture superconductivity Some of the 
typical examples with spin gaps are the spin-1 antiferro- 
magnetic Heisenberg chain ||^, the double chain spin-1/2 
Heisenberg model ||^, the spin-1/2 Heisenberg antiferro- 
magnet on a Kagome lattice Q, and Kondo spin liquid 
phase of the Kondo lattice model at half-filling . 

Recently a new system with a spin gap was found ex- 
perimentally by Taniguchi et al. for CaV409 The 
spin gap observed by magnetic susceptibility and nuclear 
magnetic resonance (NMR) measurements is A/fce — 
107K. In this paper we propose that the underlying 
lattice of the 1/5-depleted square lattice of CaV409, 
see Fig.l, favors a new type of spin disordered phase 
which may be called as plaquette resonating-valence- 
bond (PRVB) state. 

Each vanadium ion occupies a crystallographycally 
equivalent site and is surrounded by a pyramid of oxy- 
gens. First let us discuss the electronic state of this clus- 
ter (VOs)^^. In this configuration the vanadium ion is 
in V*"*" with one d-electron. Since V^"*" is surrounded by 
a pyramid of oxygen ions, the d-electron is in either dxz 
or dyz orbital. This two-fold degeneracy is lifted by a 
small Jahn- Teller distorsion whose existence is reported 
in ^ but its details are not yet clear 0. However for 
our discussion of the spin gap, the details are not impor- 
tant because has a magnetic moment of spin-1/2 for 
which single ion anisotropy is absent. 

The couplings between the spins on the vanadium ions 
are mediated by superexchange via the oxygen orbitals. 
The nearest neighbor vanadium ions share an edge of 



the square of oxgens (edge sharing), while the next near- 
est neighbor pairs share an oxgen at a corner (corner 
sharing). Superexchanges between the spins are possible 
through the hybridization with the pz orbitals of these 
oxgens. Since the number of paths for the edge sharing 
and the corner sharing is two and one, respectively, we 
expect Ji(edge sharing) = 2 J2 (corner sharing). Thus an 
appropriate model for CaV40g is the spin-1/2 Heisenberg 
model on the 1 /5-depleted square lattice 
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The magnitude of the exchange couplings may be esti- 
mated ^from the susceptibility data At high tem- 
peratures it is reasonable to assume a Curie- Weiss be- 
havior and the Weiss constant is given by k-QO = ^s{s -\- 
l)(zi Ji -I- Z2J2), where s is the spin quantum number, zi 
is the number of nearest neighbors and Z2 the number 
of next nearest neighbors. For the 1/5-depleted square 
lattice zi = Z2 = 3. From the intersect of the inverse 
susceptibility with the temperature axis it is estimated 
as 6* = 220 K Under the assumption of Ji 2J2, we 
obtain Ji/fce 200 K. 

To understand specific features of the 1/5-depleted 
square lattice let us consider a cluster of four spins on 
a plaquette. The Hamiltonian of this system is 

Hplaquctto = -^l(si + S3) • (S2 -|- S4) + J2(si • S3 -f- S2 ■ S4). 

(2) 

This Hamiltonian is readily diagonalized as is shown in 
Table I, where 513 ('S'24) is the spin quantum number of 
Si + S3 (s2 + S4) and S is the total spin quantum number. 

Table I Eigenstates of a plaquette 
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It is seen that the ground state of the plaquette is a 
singlet and the first excited state is a spin triplet with 
a spin gap of Ji for Ji/2 > J2. At Ji = 2J2 the spin 
neutral excitation crosses with the spin triplet excitation. 
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To proceed further let us first discuss a simplified 
Heisenberg model on the 1/5-depleted square lattice with 
only nearest neighbor couplings (J2 = 0). This model is 
very interesting in its own right, because the ground state 
of this model is probably disordered having no long range 
order or it is at least very close to a quantum phase tran- 
sition, in spite of the fact that the lattice is bipartite and 
completely two-dimensional. For this lattice structure, 
there are two topologically inequivalent bonds. One type 
of bonds which will be called as plaquette-bonds here- 
after form a plaquette covering of all spins, Fig. 2. We 
introduce an exchange coupling, J for this type of bonds. 
The other type of bonds, dimer-bonds, form a complete 
dimer covering of the lattice and we use J' for the ex- 
change coupling of this type. Although our final aim is 
to discuss the ground state of the model with J — J' , for 
the time being we consider J and J' as independent pa- 
rameters. In comparison, all bonds are equivalent in the 
square lattice. A dimer or plaquette covering is possible 
in this case too. However, in contrast to the above lattice 
it is not unique and it is known that the square lattice 
favors an antiferromagnetically ordered ground state. 

We begin our discussion from the limit J ^ J' . In 
this limit, as we have discussed (Table I with Ji = J 
and J2 = 0), four spins on each plaquette form a sin- 
glet ground state which has a character of resonating 
valence bond [Q. The ground-state energy of this state 
per spin is £'prvb — -~^J. In the other limit, J' 3> J, 
the global ground-state is a collection of dimer singlets, 
-E^dimcr = " |<^'- ^or J J' , another possible ground 
state is one with antiferromagnetic long-range order, 
since the lattice is bipartite. The energy of the classi- 
cal Neel state is i?Neci — ~\J' ~ \ J- These energies 
are plotted in Fig. 3. The three states discussed here 
may be considered as the simplest variational states for 
three different phases. At this level, the spin disordered 
phase, either the dimer phase or the PRVB phase, has 
a lower energy than the Neel state. Concerning the two 
singlet phases it should be mentioned that they are differ- 
ent phases although the global symmetry properties are 
the same for the two phases. It can be shown that the 
wave function of the dimer singlet state has no overlap 
with that of the PRVB state based on the fact that the 
two wave functions have different transformation proper- 
ties under the reflection with respect to the dimer-bonds: 
odd for the dimer singlet while even for the PRVB singlet. 
Thus, the two states cannot be continuously connected 
by tuning J and J' from one to the other limit. 

It is necessary to improve the estimate of the ground- 
state energies for the three phases. For the singlet phases 
we can use perturbation theories. Let us take the exam- 
ple of the PRVB state. When a dimer-bond is introduced, 
polarization processes from the singlets on both ends of 
the bond should be included. Polarization energy per 
bond may be calculated by the second order perturbation 
as — (43/576)( J'^/ J). Similar perturbation caluculation 
is also possible from the other limit of dimer singlet. The 
results are summarized as 



W = + (3) 

E^^n...^-lj'[l + \{^f] ■ (4) 

These energy are also plotted in Fig. 3. 

For the Neel ordered phase a possible improvement is 
obtained by the linear spin wave theory |^ . Extension of 
the linear spin wave theory to the present case is rather 
complicated but straightforward. The ground-state en- 
ergy per spin in this approximation is 

^^N6ci( J J') = -J^s^ + 0.325248s] , (5) 

where s = 1/2 for the present model. At J — J' the 
energy obtained by the spin wave theory, — 0.5376J, is 
very close to that of the PRVB state estimated by the 
second order perturbation, — 0.5373J. In the spin wave 
theory it is possible to calculate the reduction of magnetic 
moment by the zero point fluctuations of the magnons, 
6s = 0.288 which is nearly 50 % larger than the reduction 
for the square lattice, 5s — 0.197 and amounts to 58% 
of the magnitude of spin. The linear spin wave theory 
shows that the Neel order survives at this level but is 
on the verge of quantum phase transition. In view of 
the fact that the linear spin wave theory has a tendency 
to favor the Neel order more careful treatments are 
necessary. 

If there is a transition from the disordered phase to the 
Neel phase it would be a second order transition. In this 
case the spin gap vanishes at the transition point. There- 
fore the critical point may be estimated by examining at 
the spin gap in the disordered phases. A spin-triplet exci- 
tation in a plaquette is mobile. It can hop to a neighbor- 
ing plaquette with an effective hopping matrix element, 
J'/6, to a second neighbor plaquette with — J'^/36J and 
to a third neighbor plaquette with — J'^/216J. The po- 
larization energies for the bonds connected to the triplet 
are different form the polarization energy in the ground 
state, -(289/3456)(J'2/J). From these results the spin 
gap is calculated as 

ApRVB = J[l --(-)-—(-)'] • (6) 

Similar second order perturbation gives the spin gap for 
the dimer phase, 

^^^..^J'[l-j,-\{ j,f] ■ (7) 

Within the second order perturbation theory, the spin 
gap vanishes at {J' / J)c — 1.215 from the PRVB side and 
at {J / J')c — 0.732 from the dimer side. These points 
are shown by the dots in Fig. 3. This result suggests 
that in the narrow region between these critical points, 
antiferromagnetic long range order would exist. However, 
the spin gap remains finite at J = J' within the present 
perturbation theory. 
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An alternative way to estimate the critical points is 
a cluster mean field theory. We explain the method by 
a simple example. For the PRVB state we consider a 
cluster with four spins on a plaquettc under the influence 
of molecular fields coming from the dimcr-bonds, 

HcMF = J{si + S3) • (S2 + S4) - J'a{sl -s^2 + 4~ S4) ■ 

In the cluster mean field theory the average of a spin 
is determined by the self-consistence equation, a — (sf). 
This four spin problem can be solved analytically and the 
critical value is obtained as (J'/J)c = 3/4. When a big- 
ger cluster of 16-spins is used the critical value increases 
to (J'/ J)c = 0.8044. One can use a similar cluster mean 
field approximation for the dimer singlet. The small- 
est cluster of two spins gives (J/J')c = 1/2. The next 
one of 8-spin cluster gives (J/J')c — 0.5378. From both 
sides the critical value increases as the cluser size be- 
comes larger. Therefore we may consider (J'/,/)c = 3/4 
and {J/J')c = 1/2 as the lower limits for the critical 
points, if any. However, unfortunately, cluster size is not 
big enough to perform a reliable finite size scaling to de- 
termine the existence or absence of the Neel phase. 

All treatments discussed above suggest that the spin- 
1/2 Heisenberg model on the 1/5-depleted lattice with 
only nearest neighbor couplings has the spin disordered 
ground state in the wide region of parameter space. Only 
in a narrow region around the crossing point between the 
PRVB phase and the dimer phase, there is a possibility 
of the antiferromagnetic long range order. The recent 
quantum Monte Carlo simulations by Katoh and Imada 
| p4| suggests that there is a spin gap of A — 0.1 IJ for 
the model with J = J' consistent with the perturbation 
results. 

Let us return to the original model, Eq.(l), keeping 
the different exchanges for the dimer-bonds and the pla- 
quette bonds. It should be noted that the second nearest 
neighbor coupling is frustrating for the Neel order. It 
may be best illustrated by considering the cluster mean 
field theory. Again we consider the smallest cluster for 
the PRVB singlet, 

HCMF = J(S1 + S3) • (S2 + S4) + J"(si • S3 + S2 • S4) 

-{f -2r)a{sl~s^, + sl-sl) . (9) 

Since the eigenstates without the molecular field are 
completely determined by quantum numbers listed in 
Table I, the critical value is obtained analytically as 
( ^j'' )c = |. Thus we may conclude that the model for 
CaV409 has the quantum disordered ground state with 
safe margin. 

We can extend the second order perturbation theory 
for the spin gap of the PRVB state with the frustating 
exchange coupling, for which we obtain 

ApRVB ^J{l~hx- 2y) - ^(7x2 - lOxy + lOy^) 
3 54 

1 2x^ - ixy + Ay'^ 7 {x - yf + 

12 2-y 18 3-y ^ ' 



I + 27/2 5 (a; _ y)2 + y2 

"12 3-2y ~ 72 4-y ^ ' 

where a; = J' / J and y = J" /J. From this result it is seen 
that the gap increases as J" increases. This behavior is 
illustrated in Fig. 4 for the case J = J' . A similar 
(8)behavior is observed in the expansion around the dimer 
limit which gives, 

Adimcr - J'{1 - x-\l - ^y) - ix-2(4 - 4y + V)}. 

(11) 

It is clear that quantitative results of the perturbation 
theory is questionable at J = J' — Ji, which is seen by 
the difference between the values of the spin gap for the 
model with only nearest neighbor coupling. A = 0.205 Ji 
is obtained by the PRVB perturbation theory on one 
hand and by the Quantum Monte Carlo simulations, 
A = 0.1 IJi on the other hand |l^]. When we use the re- 
sults of the Quantum Monte Carlo simulations the mag- 
nitude of the spin gap is too small compared with the ex- 
perimental value. The present perturbation result shows 
that the spin gap increases significantly when we include 
the frustrating next nearest neighbor exchange of the or- 
der of J2 — Ji/2, which may lead to a reasonable value 
of the spin gap compared with the experiments. 

In conclusion the spin- 1/2 Heisenberg model on the 
1/5-depleted square lattice is presented as a theoretical 
model for the spin gap of CaV409. It is shown that the 
1/5 depleted square lattice is favorable for the quantum 
spin disordered phase, which may be characterized as the 
PRVB singlet. When the frustrating exchange for the 
corner sharing bonds is included it is possible to explain 
the large spin gap observed experimentally. 
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FIG. 1. A model for CaV409 of the spin-1/2 Heisenberg 
model with the nearest neighbor (solid lines) and the next 
nearest neighbor (broken lines) exchange interactions. The 
dot dashed lines show the unit cell of the 1 /5 depleted square 
lattice. 

FIG. 2. Spin 1/2 Heisenberg model on the 1/5-depleted 
square lattice with the nearest neighbor exchange interac- 
tions. Topologically there are two different type of bonds: 
plaquette-bonds J and dimer-bonds J'. 

FIG. 3. Ground state energies for the dimer singlet, the 
PRVB siglet, and the Neel ordered phase. The horizontal 
axis is X defined by J/ J' = x/{l — x). The points where the 

spin gap vanishes are denoted by the dots. Esiv is the energy 
estimated by the linear spin wave theory at J = J' 

FIG. 4. The spin gap as a function of frustrating exchange 

coupling, y — J" /J. For J" / J > 0.354 the spin gap devi- 
ates from Eq.(lO) because the minimum of the spectrum is 
different from (7r,7r). 
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